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Abstract We numerically investigate the dynamics of rotation of several close- 
in terrestrial exoplanets candidates. In our model, the rotation of the planet is 
disturbed by the torque of the central star due to the asymmetric equilibrium 
figure of the planet. We use surfaces of section and spectral analysis to explore 
numerically the rotation phase space of the systems adopting different sets of 
parameters and initial conditions close to the main spin-orbit resonant states. We 
show that, depending on some parameters of the system like the radius and mass 
of the planet, orbital eccentricity etc, the rotation can be strongly perturbed and 
a chaotic layer around the synchronous state may occupy a significant region of 
the phase space. 55 Cnc e is an example. 

Keywords Celestial Mechanics ■ Terrestrial Exoplanets • Spin-Orbit resonances 



1 Introduction 

1.1 Motivation 

Several works on physics of exoplanets assume an equilibrium state for the rotation 
of the planet. Following we give a brief list of references of different kind of studies 
on exoplanets which suppose in their models that the planet is tidally locked: 
Measurement of the rotation and oblateness of transiting giant exoplanets (e.g. 
Seager and Hui 2002; Barnes et al. 2003; Carter and Winn 2010); Generation of 
magnetic fields in exoplanets (e.g. Zuluaga and Cuartas-Restrepo 2012; see also 
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Lammer et al. 2009 and references therein) ; Atmosphere and climate of exoplanets 
(e.g. Lammer et al. 2008, Dobrovolskis 2007, Dobrovolskis 2009, Showman and 
Polvani 2011, Castan and Menon 2011); Secular orbital evolution due to tidal 
effects (e.g. Jackson et al. 2008a; Rodriguez et al. 2011); Tidal heating (e.g. Levrard 
et al. 2007, Wisdom 2008, Jackson et al. 2008b, c); Interior structure of exoplanets 
(e.g. Ragozzine and Wolf 2009, Batygin et al. 2009; Kramm et al. 2011). 

It is well know that the detected exoplanets have in general unusual orbits 
when compared to the Solar System's planets. For instance, many planets are 
much closer to their stars than Mercury is to the Sun. Some of these planets, the 
so-called close-in planets, have also eccentric orbits (see Table I). Thus, we can 
ask if the strength of the torque of the star in their rotational motion may be 
important. The main task here is to investigate this question. There exist in lit- 
erature several studies on the global dynamics of rotation of fictitious exoplanets 
which explore analytically and/or numerically the rotation phase space (e.g. Kiti- 
ashvili and Alexander 2008, Celletti and Voyatzis 2010). Our contribution here is 
to give details on the dynamics of rotation and stability of several known Earth- 
like candidates, in the vicinity of the synchronism and other spin-orbit resonant 
states. 

In this work we focus the attention in the conservative model, where we do not 
investigate the effects on the spin rate of the planet due to dissipative forces, like 
those associated to tides. Therefore we consider here only the gravitational torque 
due to the non-sphericity of the planet figure. However, our investigation can be 
useful in the interpretation of results in which tidal forces are taken into account 
in the rotational history of the planet rotation (e.g. Rodriguez et al. 2012). 



1.2 The model 

Several aspects of the physics of rotation of close-in exoplanets can be understood 
with the 1.5-degree-of- freedom dynamical system often applied in studies of phys- 
ical libration in rotation of regular satellites of the Solar system, the Moon and 
Mercury (e.g. Goldreich and Peale 1966): 

■■ 3 (B-A\ Gm . , . , > 

. y/G(m Q + m)a(l-e2) _ q(l - e 2 ) 

7 r 2 ' T 1 + ecos/' 1 > 

where: 8 is an angle measured from an inertial line, G is the gravitational constant; 
mo, m are the masses of the central star and the planet, respectively; A < B < C 
are the principal moments of inertia; a, e, / are the semi-major axis, eccentricity 
and true anomaly, respectively. It is assumed non-circular (elliptical) shape of 
the equator of the planet and null obliquity. Therefore, in the model, the rate of 
rotation of the secondary body is disturbed by the torque of the star on the non- 
spherical shape of the planet. We consider the two-body problem approximation 
for the motion of the systerrQ so that the mutual perturbations between the bodies 
belonging to multiple systems (see Section 1.4) are neglected in this paper. 



1 In this work we consider the case of prograde orbits and rotation. However, exoplanets 
in retrograde orbits with respect to the rotation axis of the star have been detected (e.g. 
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A basic assumption in the application of the model ([TJ is that the planet has a 
permanent quadrupole structure. Examples of bodies with this physical property 
in the solar system are the terrestrial planets and the regular satellites, all of which 
are, in spite of the details of their interior structures, solid-like bodies. In the case 
of exoplanets, CoRoT-7b, Kepler-lOb and 55 Cnc-e are examples of probably solid- 
like Super-Earths (Leger et al. 2009, Batalha et al. 2011, Valencia 2011, Gillon et 
al. 2012). They are transiting planets which physical parameters are being studied 
in several works. For instance, the mass of CoRoT-7b has been estimated by several 
teams obtaining different results (see Ferraz-Mello et al. 2011). 

It is still difficult to be sure about the structure of probably Earth-like exo- 
planets at the current status (see Valencia 2011, Leger et al. 2011, for the case of 
CoRoT-7b). In the case of exoplanets, the actual knowledge on planetary struc- 
ture have suffered important developments (see Baraffe et al. 2010, Valencia 2011, 
Leger et al. 2011, and references therein). In particular, Valencia et al. (2007a) 
show that, for a determined value of the mass of a rock planet in the interval 
m < WMe (Me is the Earth mass), its radius must be smaller than a determined 
value which depends on m. The limits of the radius for a given mass are the fol- 
lowing: 6600; 8600; 10,440; 11,600 and 12,200 km for 1, 2.5, 5, 7.5 and 10 M E , 
respectively. So, if a planet with m < 10M e has a radius larger than its correspond- 
ing value in the above interval, its composition can have important quantities of 
material with density smaller than the rock material, like water for instance. In 
that cases a permanent deformation is not, in principle, guaranteed. In the same 
sense Hot-Jupiters, similarly to the cases of the solar system's jovian planets, are 
probably composed mainly of gas and we also do not include them in this study 
because additional discussions are necessary in order to consider a fixed structure 
for these bodies (see Levrard et al. 2007, Rodriguez et al. 2012). 

Our main goal in this work, i.e., the application of the model ([TJ) for solid- like 
planets, has a main limitation: the total ignorance of the quantity B ^j A which 
appears in the amplitude of the perturbation of the planet rotation. In Section 2 
we propose a method to estimate the value of B (j A of Super-Earths with known 
values of radius and mass. In Section 2 we also show how to infer q as a function 
of several parameters (e.g., mass, distance, radius etc) of systems whose planets 
do not have transiting measurement and therefore do not have a determination of 
the radius of the planet. 



1.3 Methodology 

Wisdom et al. (1984) study solutions of the Equations (1 1I2[) through the com- 
putation of surfaces of sections. They numerically show the possibility of chaotic 
rotation around the synchronous and other spin-orbit states for non-spherical satel- 
lites orbiting their planets in compact and eccentric orbits. We use the Wisdom's 
methodology and apply it in the rotation dynamics of exoplanets. In Section 3 we 
show the main results of numeric integrations of Equations (1 1121) for several plan- 
ets listed in Table I. In some cases we compare the results obtained with surfaces 

Hebrard et al. 2011, and references therein). Retrograde rotation is also a possible equilibrium 
state (Correia et al. 2008). The study of spin-orbit resonances in retrograde motion should be 
interesting. 
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of section with dynamical maps based in spectral analysis of the solutions (e.g. 
Michtchenko and Ferraz-Mello 2001). 

Our simulations cover wide ranges of the free parameters including moments of 
inertia ratio and orbital eccentricity. We also give the range of some other param- 
eters (like planetary radius, mass) which would correspond to solid-like planets 
in view of structure models of exoplanets (e.g. Valencia 2007a, b; 2011). Thus, we 
discuss the dynamics of planetary rotation with our model in large sets of ini- 
tial conditions and appropriate parameters corresponding to solid-like planets. We 
consider planets with orbital period P < 33 days and mass m < 15M^. Table 
I shows the data of several Super-Earths considered here. Inspection of Table I 
helps us to see some common properties of all planets which we study here. For 
instance, there are several planets with published value of non-circular orbit, even 
with very small values of orbital period P < 7days. As we see in Appendix A (see 
below), the existence of non-circular orbits for close-in planets is in conflict with 
the classical results of orbital circularization of close-in orbits due to tidal effects. 



1.4 Orbital eccentricity 

As we will see in Section 3, the planets with more rich rotational dynamics are those 
which move in eccentric orbits. In fact, the planet's orbital eccentricity plays an 
important role since in the case of non-circular orbit the amplitude and argument of 
the sine in (JXJ) suffer temporal variation. However, the existence of eccentric close- 
in orbits in single-planet systems cannot be guaranteed. In fact, an initial eccentric 
orbit would be tidally damped, evolving to circularization in timescales depending 
on physical parameters of the interacting pair and also on the initial values of the 
elements (see Dobbs-Dixon et al. 2004, Ferraz-Mello et al. 2008, Rodriguez and 
Ferraz-Mello 2010, Rodriguez 2010). In Appendix A we calculate, using averaged 
equations, the circularization timescale of the Super-Earth CoRoT-7b, adopting 
different dissipative parameters. Our results, in agreement with previous ones, 
show that tidal evolution of systems with close-in companions with orbital period 
less than one week can be the responsible for the circularization of the orbits. 

Thus, probably the values of the eccentricities of some close-in single-planets 
listed in Table I are not well determined by the current methods of orbital fitting, 
and/or some additional mechanism can explain the non-null values. For instance, 
under the hypothesis of an additional undetected resonant companion, a large value 
of eccentricity of a single planet can be in fact close to zero, a result that shows a 
degeneracy of the problem of orbital determination of the system (Anglada-Escude 
et al. 2010). Tadeu dos Santos et al. (2012), in a paper written with the main 
propose of study the dynamics and the determination of the orbital elements of Gl 
581g, show degeneracy in the determination of the eccentricity of all members of 
the Gl 581 system. In this case, degeneracy means that depending on the planetary 
architecture adopted in the model (e.g. number of planets), different solutions are 
possible with the same statistical significance for different values of eccentricity. 

The existence of exterior companions can also work as a dynamical mechanism 
which may explain the non-null eccentricity values of close-in planets of multi- 
planetary systems. Spiegel et al. (2010) evaluate that terrestrial planets at 1 AU 
can be highly perturbed by an outer giant planet, which may excite large variations 
of eccentricity of the inner planet. At such distances, however, tidal effects may be 
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negligible (Appendix A), and eccentricity may be large. However, previous works 
have show that, the combination of secular interaction and tides raised on the 
inner planet by the star result in the circularization of both orbits (see Mardling 
2007; Rodriguez et al. 2011). Indeed, depending on the planet mass ratio, the 
perturbation of the companion acts to increase the circularization timescale of 
the inner orbit. Hence, the explanation for the non-null eccentricities of close-in 
planets remains under discussion (see Matsumura et al. 2008, Correia et al. 2010, 
Correia et al. 2012). 

Therefore, notwithstanding all the above discussions, we consider in this work 
the orbital eccentricity as an open parameter in our calculations. The reported 
values, listed in Table I with the error bars, will also be taken into account for 
sake of reference. 



Table I: Data of planets with orbital period P < 33days (except Kepler-llf ) and 
mass m < 15Me, where Me is the Earth mass. The list is sorted in increasing value 
of the orbital period. The data have been taken from http://exoplanet.eu/index.ph 
except in other cases indicated. The other units are the solar mass (Mst/jv) and 
the Astronomical Unit (AU). a : Ferraz-Mello et al. (2011); b : Planet present in 
multiple system; c : See also Section 3.2.1. d : Value taken from Charpinet el al. 
(2011), differing from the value given in http://exoplanet.eu/index.php. 



Planet m: Mass mo: Star's P: Orbital a: semi-major e: 

(Me) mass (Mgupf) Period (day) axis (AU) eccentricity 



KOI-55b 


0.445 d 


0.496 


0.2401 


0.006 


? fc 


KOI-961c 


1.907 


0.13 


0.4533 


0.006 


9 b 


KOI-55c 


0.667 


0.496 


0.3429 


0.0076 


? b 


55 Cnc e 


8.58 


0.905 


0.7365 


0.0156 


< 0.06 6 


Kepler- 10b 


4.55 


0.895 


0.8375 


0.0168 


0» 


CoRoT-7b 


8.5 a 


0.93 


0.8536 


0.0172 


b 


KOI-961b 


2.8604 


0.13 


1.2138 


0.0116 


? 6 


GJ 1214b 


6.357 


0.153 


1.58 


0.014 


< 0.27 


Kepler-9d 


6.99 


1.0 


1.5928 


0.0273 


? 6 


KOI-961d 


0.9535 


0.13 


1.8516 


0.0154 


76 


GJ 876d 


6.67 


0.334 


1.9378 


0.02 


0.207 ± 0.055 6 > c 


GJ 3634b 


7.0 


0.45 


2.6456 


0.0287 


0.08 ± 0.057 


Kepler-21b 


< 10.488 


1.34 


2.7857 


0.0425 





Kepler-18b 


6.897 


0.972 


3.5047 


0.0447 


0» 


Kepler-20b 


8.58141 


0.912 


3.6961 


0.0454 


< 0.320 6 


CoRoT-7c 


8.39 


0.93 


3.698 


0.046 


0» 


HD 156668b 


4.16 


0.772 


4.646 


0.05 





GJ 674b 


11.76 


0.35 


4.6938 


0.039 


0.2 ±0.02 


HD 7924b 


9.217 


0.832 


5.3978 


0.057 


0.17 ± 0.16 


HD 45184b 


12.7132 




5.8872 


0.0638 


0.3 


Kepler-20e 


< 3.083 


0.912 


6.0985 


0.0507 


■fb 


GJ 176b 


8.42 


0.49 


8.7836 


0.066 





HD 97658b 


6.356 


0.85 


9.4957 


0.0797 


0.13 ±0.06 


HD 125595b 


14.30 


0.76 


9.67 


0.078 





Kepler-llb 


4.3 


0.95 


10.3037 


0.091 


() b 


Kepler- 11c 


13.507 


0.95 


13.025 


0.106 


o b 


Kepler-20f 


14.3 


0.912 


19.577 


0.11 


fb 


Kepler-lld 


6.099 


0.95 


22.687 


0.159 


>' 


Kepler- lie 


8.40 


0.95 


31.9959 


0.194 


o b 


Kepler-llf 


2.3 


0.95 


46.6888 


0.25 
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2 The Prolateness 

The first quantity which we must have in hand in order to study the system 1111211 is the 
ratio B q A ■ Here we propose a method to estimate the value of B q A as a function of several 
parameters of the system, like the mass of the star, the mass of the planet, the semi-major 
axis, and the radius of the planet. 

Consider a homogeneous non-viscou^ rotating planet orbiting a star. The equilibrium 
figure of a rotating body (with axes a, b and c) under the action of the tide generating force 
of the central star, lying in the equatorial body's plane, is given by a Roche ellipsoid in which 
a > b > c (see Chandrasekhar 1969). The smaller axis is perpendicular to the largest axis, 
which is pointed in the direction of the central body in the case of synchronous rotation. In this 
case, we have the following expression for the deformation of the equator, or the prolateness, 
hereafter denoted by e: 



_ a-b _ 15 wo (R\ 3 ( 5(mo+m) /ft\ 2 ( R\ 3 \ 
£_ a ~~ 4 m \r) [ 2m \n) \r) J " 

(e.g. Ferraz-Mello et al. 2008). 

When the rotation is neglected, the attained equilibrium figure is a Jeans spheroid with 
b = c. In this case, the value of e is given by equation above without the second term. A 
classical result shows that for a free synchronous rotating body the polar deformation or 
oblateness (that is, 1 — r) is three times smaller than the prolateness due to the companion 
(when rotation is neglected, see Danby 1988, p. 121). The restriction b = c enable us to write 
e in a convenient form. Indeed, the momenta of inertia of the spheroid can be written as 



A=-m(b 2 +c 2 ), B=im(a 2 + c 2 ), C = -m(n 2 + b 2 ). (4) 
5 5 5 

Hence, using the definition of e and imposing b = c into above equations we arrive to 

A 2 1 



1-6- -6^. (5) 



C (l-e)- 2 + l " 2 
Therefore, at first order in e: 

B — A 

~ e. 6 

C y ' 

Contrary to the case of some bodies of the Solar systerr0, the values of the momenta of 
inertia of exoplanets are completely unknown quantities. However, we can take an idea of the 
magnitude of the prolateness by using equations (0 and J5J . 



2 In this hypothesis we are neglecting the dynamic tidal torque raised by the central star 
which also affect the planet rotation. However, several planets which we study here are so 
close to their star that tidal effects may be important, mainly in the case of eccentric orbits. 
We devoted to Appendix B the illustration of one case in which the tidal torque is taken into 
account and applied to CoRoT-7b super-Earth planet. 

3 In order to test the validity of Equation ||6} we estimate e for some bodies of the Solar 
System. For Europa we estimate e ~ 1.9 X 10~ 3 , a value for which the order of magnitude 
agrees with current estimative of ~ 1.5 X 10 — 3 . The latter value has been obtained from the 
numerical values of the gravitational coefficient C22 and mean moment of inertia £ = c „„ 
given in Van Hoolst et al. (2008) through the relation e = 4C22/£- In the case of Titan, 
we obtain C22 = 1-4 X 10 — 5 , in agreement with one reported value of about 1 X 10~ 5 (less 
et al. 2010). For Io, we find C22 = 3.4 X 10~ 4 , in concordance with Schubert et al. (2004), 
who reported C22 = 5.6 X 10 — 4 . For Mercury we obtain e ~ 10~~ 6 , while its current order of 
magnitude is ~ 10 — 4 . Mercury has an eccentric orbit (0.205), but our estimative of e does not 
change considering Mercury at the pericenter. Probably the dependence of e on the eccentricity 
and distance be more complex in the case of more distant and eccentric planets. 
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The dependence on (^j in {[3} is in fact weak. For instance, Table II shows the values 
of e considering Q = n for several transiting Super-Earths which has a determination of the 
radius. In the case of 55 Cnc e, e~2.8x 10~ 2 . For £2 = 2n, e ~ 2.97 X 1(T 2 . If we neglect 
the rotation term and consider Q = 0, we have e ~ 2.75 X 10 — 2 . These values are smaller than 
when the uncertain in radius of planet is considered in the calculation of t. For instance, in 
the case of 55 Cnc e, Gillon et al. (2012) give R = 2.17±0.10 Re, which results the maximum 
values: e ~ 3.21 x 1CT 2 (Q = n), e ~ 3.43 x 10" 2 {Q = 2n), e ~ 3.14 x 10~ 2 (17 = 0). 

All other planets listed in Table I do not have an estimative for R. In these cases we can 
estimate the order of magnitude of e by analyzing its dependence with different parameters, as 
we illustrate below. Figures l(a-e) show the level curves of the function e = e(r, m) calculated 
with Equation (0 considering five fixed values of the radius of a fictitious planet of a star with 
mass mo = ImsijN, where msuN is the solar mass. The contours given in gray color corre- 
spond to "small" value of e < 0.01, while the colored levels correspond to "large" prolateness 
where e > 0.01. We arbitrarily define here the value e = 0.01 as the limit between large and 
small prolateness of Super-Earth^|. Figure 1(a) shows that exoplanets similar to the Earth 
(i.e., with mass and radius similar to the Earth's values) can have large equatorial prolateness 
for all r < 0.018AU (note the colored level curves near to the bottom left of the panel). 

A zoom close to the origin of the Figure 1(a) is shown in Figure 1(f). We can note that for 
small values of (r, m) in the figure we have e ~ 1 (blue levels J^. Thus Figure 1(f) suggests that 
a planet with r 0.018 can suffer important effects of the torque on its rotation. The recently 
discovered members of the planetary system of the stars KOI-961 and KOI-55 are examples 
of such kind of planets orbiting sub solar mass stars (see Table I). As they are transiting 
planets we can calculate their e. Table II shows that the prolateness of KOI-55b, KOI-55c and 
KOI-961c are large, i.e., e > 0.01 while in the case of KOI-961b,d, e < 0.01. In Section 3.1.2 
we will study the rotational dynamics of these planets in details. 

Finally we note that in the present study we have assumed constant B ^ A value. However, 

depending on the adjustment to hydrostatic equilibrium B q A may be itself time-dependent. 
The time-scale for adjustment to the completely relaxed state depends on material parameters, 
e.g., rigidity and viscosity which are highly temperature dependent. Furthermore, it will depend 
on the spin-orbit coupling the planet is trapped in. For the application to terrestrial planets 
given here the constant, i.e. frozen, B ^ A configuration is justified because adjustment of solid 
rocky planets to equilibrium usually occurs on geological time-scales. Therefore, the states 
described here do not necessarily represent hydrostatic configurations. For Jupiter-like planets 
the simplification of constant — ^A- values may not be justified because of the fluid behavior 
of the gas giants even on the short time-scales of their orbital periods. Adjustment to the 
hydrostatic state in a specific resonance will be relatively fast. In this case the time-dependence 
of the gas giants B ~q A must be taken into account in a coupled set of differential equations 
describing the spin, the orbit and the hydrostatic gravity field of the planet. 



4 As we will see in details in Section 3, for e > 0.01 the rotation may be highly disturbed 
by the torque of the central star in the case of eccentric orbits. However, this division between 
"small" and "large" prolateness is arbitrary and will be used in this work only for didactic 
purposes. 

5 We remark that very large values of e close to the unit are not valid within the first-order 
approximation in e. 
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Table II: Numerical value of B q A — £ (Equation (1)) for several Super-Earths 
with known radius, e has been calculated with Equation fl3j where we consider 
Q = n. The values of radius, in units of the Earth radius (IRe = 6,378 km), have 
been taken from http://exoplanet.eu/index.php, a : Charpinet el al. (2011); b : The 
current estimative in http://exoplanet.eu/index.php is 2.84±0.2 Re; c This value 
differs from those given in Version 1 (vl) (e = 0.117) because we use the mass of 
the planet given in Charpinet el al. (2011), which is ten times smaller than the 
value given in http://exoplanet.eu/index.php. 
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Fig. 1 (a-e) Level curves of e in the plane (r X m) (star-planet distance (r) and the mass of 
the planet (m)), for different values of the radius of a fictitious planet orbiting a star with 
1 Msun- (f) Zoom of Figure 1(a) close to the origin of the r-axis. The blue contour shown 
in (f) shows prohibitive levels with e ~ 1. The white portion of the figure close to the origin 
corresponds to the levels curves which are not physically possible to exist since they correspond 
to e > 1. The units are the Earth radius (Re), the Earth mass (Mg), the solar mass (Msun) 
and the Astronomical Unit (AU). 
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3 Surfaces of Section and the Dynamics of Rotation 

In this section we show surfaces of section of the numerically computed solution^] of the model 
iffll for several Super-Earths listed in Table I. Following Wisdom et al. (1984), we investigate 
the dynamics of rotation in the plane (8, 8/n) at / = 0, where / is the true anomaly and 8/n is 
the ratio between the angular rotational velocity and the mean motion. Thus, in the surfaces 
of section, we plot the values of (8,8 /n) at the moment of the passage of the planet through 
the pericenter in the case of eccentric orbit. 

In order to present surfaces of sections, we consider four main types of Super-Earths: 
planets with circular or non-circular orbit, and those which are single planets or are present in 
multiple systems. We divide the presentation of the results as follows. In Section 3.1 we review 
basic topics of spin-orbit resonances and study examples of close-in planets (single or not) 
with large equatorial prolateness and circular orbits. In Section 3.2 we study the dynamics of 
rotation of close-in planets present in multiple systems considering eccentric orbits with orbital 
period P < 7.5 dayfl In Section 3.3 we analyze the cases of single planets in both, circular 
and eccentric orbits. In Section 3.4 we apply the Chirikov overlapping criterium to explore 
chaos in some systems around the domains of the 1:1 and 3:2 spin-orbit resonances. 



3.1 Close-in Earth-like planets with circular orbits 

3.1.1 CoRoT-lb, Kepler-lOb and the basic topics of spin- orbit dynamics 

Figure 2 shows surfaces of sections for CoRoT-7b and Kepler-lOb, two Super-Earths very close 
to their host stars (a ~ 0.017 AU). Figure 2(a) corresponds to the case of CoRoT-7b. Since 
its orbit is circular, the system JTJ), J2]l is integrable in this case and therefore we say that 
the dynamics of rotation of CoRoT-7b is regular, i.e., no chaotic motion is admissible in the 
domain of the model. We indicate in Figure 2(a) the domain of the synchronous motion located 
around the equilibrium point at 9 = 0, 8 = n. The closed curves around the synchronous state 
show the amplitude of the physical libration of the major axis of the planel[j relative to the 
star, which motion is given quantitatively by the oscillation of the synodic angle tp = f — 8 
around zero. In the section 8 follows ifi since we put / = 0. 

In Figure 2(b) we consider a fictitious non-null value for the orbital eccentricity, e = 0.01. 
In this case, the domain of other spin-orbit resonances may appear in the sections, like the 
3:2 resonance for instance. In fact, it can be shown that the domains of low-order spin-orbit 
resonances exist in the phase space only in the case of e ^ 0. This can be understood as follows. 
The average, with respect to the mean-anomaly, of the torque of the star on the rotation of 
planets in circular orbits, is null for all low-order spin-orbit regimes except the synchronous 
one (Goldreich and Peale 1966). Thus, if the torque is null in a determined resonance, libration 
of the corresponding synodic angle is not possible. Therefore, in the case of circular orbit this 
is a physical explanation for the absence of resonant islands in the surfaces of section except 
the synchronous one. 

The domain of the synchronous motion is always present in the phase space even for e = 0. 
For e = the average torque is never null except if B — A = or at the equilibrium point 
of the exact synchronism. For circular orbits the fixed point of the synchronous motion is 
located at 8 = n and 8 = 0. For e ^ 0, the fixed point corresponding to the equilibrium 
of the system is not located at 8 = n because the existence of a forced component (re), and 
it is shifted compared to the case of null eccentricity. This can be seen by inspection of the 



6 We use the code of Everhart (1986). 

7 In the Section 4 we list the results for some planets with 7.5 < P < 33 days. 

8 Following classical astrodynamics texts (e.g. Danby 1988), we use in this work the term 
physical librations for those motion associated to the non-spherical shape of the secondary 
body. The often called optical librations do not require an asymmetry of the secondary since 
they are not associated to a disturbed rotator. For instance, the optical libration in longitude 
exists due to differences between a (constant) velocity of rotation and the instantaneous orbital 
velocity in an eccentric orbit. In the disturbed (non-spherical) case the optical libration is 
superposed to the forced component. 



Dynamics of Rotation of Super-Earths 



11 



CoRoT-7b: (B-A)/C~0.00992 



e=0.01 




Kepler-10b: (B-A)/C~0.0115 



e=0.01 




Fig. 2 Surfaces of Section in the plane (8/nx8) of numerical solutions of JTJ . The sections have 
been done at each orbital revolution of the planets CoRoT-7b (Figures 2(a,b)) and Kepler-lOb 
(Figures 2(c,d)). Each orbit has been integrated for 1140 orbital periods of the planet. The 
numerical values of B q A has been calculated with Equation 10 and are given at the top of the 
panels (see also Table II). In the plots given at left, the orbital eccentricities correspond to the 
real values (i.e., the reported ones). In the plots given at right, the values of the eccentricities 
are fictitious and have been chosen in order to evaluate the domain of the 3:2 spin-orbit 
resonance in the case of quasi-circular orbit. 
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KOI-55c: (B-A)/C~0.117 




Fig. 3 Surfaces of Section in the plane (0/n X 9) of numerical solutions of ||T)|. The sections 
have been done at each orbital revolution of the planet and KOI-55c. Different fictitious values 
of the orbital eccentricity (indicated in each plot) have been considered. Each orbit has been 
integrated for 1140 orbital periods of the planet. The numerical values of B q A has been 
calculated with Equation (0 and are given at the top of the panels (see also Table II). 
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fixed point of the 1:1 resonance in the surfaces of sections, mainly in the case of high values 
of eccentricity (e.g. Figures 6(d) and 7(c,d)). At first order in eccentricity, it can be shown 

that: = n2 "^a e ' where uJq = 3n 2 B ^ A for the 1:1 resonance (Murray and Dermott 1999, 

pages 215-16). In the surface of section, it can be shown that the equilibrium point of the 1:1 

resonance is located at — ~ 1 — re. 

For small eccentricity it can be shown that physical libration of ij) = Vi:i = / — includes 
two component: the often called forced term, with frequency given by the mean-motion n, 
and long-term, free libration, with frequency wq. The general expression of ujq is given by 
lOq = 3n 2 \H(p, e)| B q A , where: H(p,e) are the coefficient of the mean torque at the resonance 
of order p; p is a half-integer which defines the resonance such that p = 1 for 1:1 resonance, 
p = 2 for 2:1, p = +3/2 for 3:2 etc (see Appendix B and Appendix in Rodriguez et al. 2012JJ 
For 55 Cnc e, close to the 1:1 fixed point, the long-term period is about 2.54 days, while the 
forced component libratcs at each ~ 0.736 days. In the case of non-circular orbits, chaotic 
rotation is physically possible mainly close to the scparatrices of the spin-orbit resonances. 
This regime of motion is possible since the physical libration can become unstable for large 
amplitude. In the example given in Figure 2(b), however, the values of eccentricity and B ~^ A 
are very small, and the chaotic layer is very thin. 

Figures 2(c,d) show the case of Kepler-10 b. We can note the similarities with the system 
of CoRoT-7b discussed above. 



3.1.2 KOI-55b, KOI-55c 

In spite of regular motion analyzed in Figure 2, we can ask if a fictitious planet closer to its 
star than CoRoT-7b and Kepler 10-b would have similar properties of the rotation phase space 
even in the case of near-circular orbits. As discussed in Section 2, depending on the star-planet 
distance, a planet can admit very large value of e > 0.01. We give two examples of this kind 
of Earth-like candidates recently discovered: KOI-55b,c, located at 0.0060 AU (planet b) and 
0.0076 AU (planet c) from the star, respectively (Table I). The values estimated for B @ are 
~ 0.26 (KOI-55b) and ~ 0.117 (KOI-55c) (Table II). 

Figure 3 shows surfaces of section for KOI-55c. The numerical values of the orbital ec- 
centricity is unknown (Charpinet el al. 2011). Due to their proximities to the star, it can be 
expected that it eccentricity is close to zero (see Appendix A). However, we consider in Figure 
3 three different quasi-circular orbit. The domain of the 1:1 resonance is very large. The phase 
space shows the presence of other spin-orbit resonances, like the 3:2 resonance. The fixed points 
of all resonances but the synchronous one are slightly forced from the exact commensurability. 
For instance, in Figure 4(a) we show the physical libration of V3/2 = 2f — 38 adopting an 

initial condition within the 3:2 resonance indicated in Figure 3(b), where 9q = 1.685n and not 
80 = 1.5n as would be expected. 

The phase space also shows secondary resonances within the synchronous island. Sec- 
ondary resonances occur when the period of libration of the physical libration around the 
synchronous state is commensurably with the mean-motion of the planet. Figure 3(a) shows 
the 2:1 secondary resonances. 

Figure 3(c) shows that setting e = 0.01, both, the chaotic layer around the separatrix of 
the 1:1 resonance, and the domain of 2:1 secondary resonance, have large domains in the phase 
space. An orbit close to the 2:1 secondary resonance is shown in Figure 4(b). Note that, while 
the true anomaly and the rotation angle evolve at approximate same rate (period=0.34 days), 
the angle i/> = / — 8 associated to the 1:1 resonance librates with twice of the period of motion. 

The dynamics of rotation of KOI-55b is very similar to the case discussed above. 



9 Note that uiq is a linear function of n with linear coefficient depending on the prolateness 
and H(p,e). In the case of 1:1 resonance, at first order in eccentricity the relation between 
cjq and n is invariant to the value of eccentricity since H(l,e) = 1. This does not occur for 
large eccentricity for other low-order resonances like the 3:2 one, where H(3/2, e) = at first 
order. 
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Fig. 4 (a) Time evolution of V'3/2 = 2/ — 3£>. The initial condition corresponds to the small 
island within the 3:2 resonance indicated in Figure 3(b): 6>o = 5°, 80 = 1.685n, /o = 0°, 
where n is the mean motion. The long-term oscillation is the free libration while the short- 
term libration includes the forced one (see also Section 3.1.1). (b) Time evolution of different 
quantities indicated in the figure. The initial condition corresponds to the larger island of 
secondary resonance shown in Figure 3(c): 6q = 49°, 6q = n i fo = 0°. 
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KOI-961c: (B-A)/C~0.018 KOI-961b: (B-A)/C~0.0013 
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Fig. 5 Surfaces of Section in the plane (6/n X 0) of numerical solutions of JTJ. The sections 
have been done at each orbital revolution of the planets KOI-961c (a), KOI-961b (b), CoRoT- 
7c (c,d). Each orbit has been integrated for 1140 orbital periods of the planet. In (c,d), the 
numerical values of B q A has been calculated with Equation 10 supposing that CoRoT-7c has 
the same radius of CoRoT-7b. Different values of the eccentricity have been used in the case 
of CoRoT-7c. 
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3.1.3 The system of KOI-961 

The star KOI-961 hosts three transiting sub-Earth-sizcd planets (Muirhead et al. 2012). We 
estimate the B @ A of these planets (Table II). KOI-961c has the larger e ~ 0.018; its dynamics 
of rotation (Figure 5(a)) is very similar to the previous cases (e.g. Figure 2). The other two 
planets have relatively small values of e <C 0.01 and the domain of the 1:1 resonance is very 
thin in their rotation phase space, like Figure 5(b) illustrates. 

3.1.4 CoRoT-7c and the planets of the stars Kepler-18, 20 and 21 

Earth-like candidate planets in circular orbits close to the star with semi-major axes larger 
than the previous cases have regular dynamics of rotation since their prolateness are small and 
the torque of star is negligible. The domain of the phase space of the synchronous resonance is 
very thin, like we show in Figure 5(c) in the case of CoRoT-7c, where we have supposed that 
its radius is equal to the radius of CoRoT-7b in order to estimate B q A ~ e = 0.0005. 

In Figure 5(d) we put a fictitious value of orbital eccentricity of CoRoT-7c equal 0.1. A 
characteristic of the phase space in these cases is that the domains of the 3:2 and 1:1 resonance 
have similar width in the phase space. 

Supposing circular orbit, analyses similar to that given in Figure 5(c) hold to Kepler- 
21b, Kcplcr-18b and Kcplcr-20b, close-in planets with orbital parameters with same order of 
magnitude of CoRoT-7c's (see Table I). However, Kepler-21b, Kcplcr-18b and Kepler-20b are 
transiting planets and we can give an estimative of their values of B q A , which are listed in 

Table II. Table II include also the B q A of Kepler-20e and Kepler-20f, two planets present in 
the Kepler-20 system which have masses and orbital periods whose values are given within the 
range considered in this work (see Table I). 

3.1.5 Planetary System of Kepler- 11 

The value of e for the close-in planets of the star Kepler-11 are listed in Table II. Their values 
are small and the dynamics of rotation the planets are regular in our model. 



3.2 Planets with eccentric orbits with P < 7.5 days present in multi-planetary 
systems 

3.2.1 GJ 876d and Kepler-9d 

Kepler-9d and GJ 876d are probably arc small planets belonging to multiple systems where the 
outer members are hot Jupiters and whose orbits are close to the 2/1 mean-motion resonances. 
GJ 876d and Kepler-9d arc two planets with similar orbital period and mass. They orbit around 
stars with very different values of mass since GJ 876 is a M4 V star with mass ~ 0.334Msun 
(Corrcia ct al. 2010), while Kepler-9 is a solar- type star with mass mo ~ IMsun (Torres et 
al. 2011). Kepler-9d is a planet whose radius is determined but the orbital eccentricity is not. 
On the other hand the eccentricity of GJ 876d is known (while its radius is not). Therefore 
we consider the value of e as a free parameter in the case of Gliese 876d, and the orbital 
eccentricity in the case of Kepler-9d. 

If GJ 876d is a terrestrial-like, solid planet, Valencia et al. (2007b) show that its radius 
cannot be larger than ~ 12, 000 km ~ 1.88Re- Adopting this value of radius we estimate the 
prolateness of this planet: B q A ~ 0.0036. GJ 876d has one of the largest reported orbital 
eccentricities among all the close-in Earth-like candidates. In order to study its rotation with 
our model we use the values of the eccentricity given in Correia et al. (2010), e = 0.124, and 
the value e = 0.207 given in Rivera et al. (2010). The results are shown in Figures 6(a,b). 
Within the above scenario we see by inspection of the surfaces of section that the rotation of 
Gliese 876d is probably regular within the domain of our model. It is interesting to note in 
this figure that the domains of 2:1 and 3:2 resonances have similar widths in the phase space. 
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Gliese 876d: (B-A)/C~0.0036 
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Fig. 6 Surfaces of Section in the plane (8/nx8) of numerical solutions of JT). The sections have 
been done at each orbital revolution of the planet Gliese 876d (figures (a-d)) and Kepler-9d 
(e,f). Each orbit has been integrated for 1140 orbital periods of the planet. B q A of Kepler-9d 

is given in Table II; in the case of Gliese 876 d the value of B q A is discussed in Section 3.2.1. 
In (a,b) the orbital eccentricities correspond to two published values (see references in Section 
3.2.1), while in (c,d) the eccentricities utilized in the simulations are fictitious. 
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The radius of GJ 876d is unknown, so we also study its dynamics with a larger value than 
that given above: R ~ 3Re — 19,000 km. From the plots similar to those shown in Figure 
1 (Section 2), we can see that for this radius e ~ 0.01. Figures 6(c,d) show that the rotation 
can be more complicated in this case. For instance, we can see the overlap of the separatrix 
of the 3:2 and 1:1 resonances, engulfing the 5:4 resonance, when e = 0.207 (Figure 6(d)). All 
resonances, including the island of 1:2 resonance, have larger and similar domains in the phase 
space in the case of large e (compare Figures 9(a,b) and Figure 6(c,d)). 

The eccentricity of Kepler-9d is unknown, but in analogy to the case of GJ 876d, we study 
the dynamics of its rotation in the case of eccentric orbit. We investigate the dynamics of 
Kepler-9d for two values of eccentricity: e = 0.1 (Figure 6(e)) and e = 0.2 (Figure 6(f)). The 
rotation of Kepler-9d is regular in both cases. 

The similarities between Figures 6(a,b) and Figures 6(e,f) are evident. The value of e 
estimated for Kepler-9d and GJ 876d which have been utilized in these figures are similar: 
e ~ 0.003; we recall that the e of GJ 876d was calculated supposing R ~ 1.88Re given in 
the theory of Valencia et al. (2007b). Adopting the mass-radius relation of Valencia et al. 
(2007a), Kepler-9d is probably a solid-like body and our results on its dynamics of rotation are 
confident in this case. The case where we consider e = 0.015 for GJ 876d must be understood 
as a first view of a more complicated scenario since in this case the planet is not a solid-like 
body (Valencia 2007b). 



3.2.2 55 Cnc e 

55 Cnc e is a low-mass planet which has the smaller orbital period among all planets orbiting 
a solar-like star (G8-K0). With the recent estimative of its radius, it is probably that this 
body is an Earth-like planet with important rock-component in its structure (e.g. Gillon et al. 
2012 'P^l . a result which motivate us to study some aspects of its rotation with our model. 

55 Cnc e has the second larger value of e which we have calculated for all Super-Earths 
considered in this work (Table II). Figure 7 shows four plots for this planet where we adopt 
different values for the orbital eccentricity. The dynamics of rotation is essentially regular in the 
immediate vicinity of the resonances for e < 0.01. For larger values, like the case corresponding 
to the larger published value of the eccentricity (e = 0.057), is shown in Figure 7(c). In this 
case there is a large chaotic region around the separatrix of the 1:1 resonance. For larger values 
of eccentricity, a large region of chaos located between the 1:1 and 3:2 appears in the phase 
space (Figure 7(d)). Note also the presence of 1:4 secondary resonance inside the synchronous 



In order to confirm and to explore in more details the properties of the phase space given 
by the surfaces of section we compute, in Figure 7(e), a dynamical map in a grid of initial 
conditions corresponding to range of the axis given in Figure 7(c) . In the dynamical maps we 
plot the spectral number, N, defined by the number of peaks given in the numerical Fourier 
spectrum of some variable of the problem which are greater than a fraction of the highest 
amplitude (e.g. Callegari and Yokoyama 2012). We adopt reference amplitude of 10% of the 
largest peak present in the spectrum of 8. N is represented on each point of the grid by a 
color given in a gray scale which value depends of N in the following way: i) white points 
corresponding to N = 1; ii) gray points where 1 < N < N c , with N c being a (arbitrarily) 
large value adopted for N; iii) black points with N > N c . We adopt N c = 100 in Figure 7(e). 
Distinct regions are different stability regions of the phase space: in white and gray regions 
the motion is regular, where the spectra of the solutions contains a small number of separated 
and countable peaks. In darker regions, where the spectrum contains a large number of peaks, 
the motion can be chaotic, and in general the peaks are not well separated in the spectrum. 
Comparison and inspection of Figures 7(c,e) show us that the dynamical map of frequencies 



10 According to Valencia (2007a), supposing a solid-like structure for 55 Cnc e, its radius 
cannot be larger than ~ 1.8Re- Since the estimated radius of 55 Cnc e is ~ I.YIRe it cannot 
be considered a pure solid-like planet, admitting a non-negligible layer of other components 
like water in its composition. Thus, it can be interesting to improve our model to consider the 
perturbations in the rotation of planet adopting different layers (e.g. van Hoolst et al. 2008). 

11 Wisdom (2004) shows that, in the case of regular satellites of the Solar System, a capture 
in such rotational state may enhance the tidal heating in the satellite interior by several orders 
of magnitude. The same would occur for close-in Super-Earths. 
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agrees very well with the main dynamical regions of the system obtained with the surfaces of 
section (i.e., spin-orbit resonances, chaotic layers, secondary resonances). 

We conclude that a non-null value of the eccentricity of 55 Cnc e, and its proximity to the 
star, lead a strong periodic perturbation on its rotation. Thus, according to our model, in a 
evolutionary scenario where some dissipative mechanism is able to reduce the primordial spin 
of the planet to the synchronous resonance, 55 Cnc e would have crossed a chaotic layer before 
the synchronism be attained, similar to the history of the rotation dynamics of the irregularly 
shaped, close-in, satellites of the solar system (Wisdom et al. 1984, Wisdom 1987, Wisdom 
2004). This conjecture has been confirmed by Rodriguez et al. (2012) in their scenario of tidal 
evolution of the system of 55 Cnc. 



3.3 Systems with a single planet 

We consider in this separated section some aspects of the dynamics of rotation of systems with 
a single planet. Consider the cases of planets with circular and eccentric orbits separately. 

The systems with single planets in circular orbit which we study are: HD 156668b (orbital 
period P = 4.646days, m = 4.16M B ), GJ 176b (P = 8.78days, m = 8A2M E ) and HD 
125595b (P = 9.67days, m = 14.3Mb). The radius of all of them are unknown quantities, 
but we can use Figures 4(g,h,i,k) to infer that their prolateness are negligible due to their 
relatively large mass and planet-star distance (a > 0.05AU). Therefore their rotations suffer 
very small perturbation of their respective stars. In the case of HD 156668b the prolateness is 
slightly larger than the other cases due to its relatively smaller mass and distance to the star, 
but anyway its rotation is regular since we assume circular orbits. 

There are several examples of isolated Super-Earths with reported eccentric orbits: GJ 
3634b (P = 2.64days, e = 0.08, m = 7.0M E ), GJ 674b (P = 4.7days, e = 0.1, m = 11.76M B ), 
HD 7924b (P = 5.4days, e = 0.17, m = 9.27M B ), HD 45184b (P = 5.88days, e = 0.3, 
m = 12.7Mb), HD 97658b (P = 9.5days, e = 0.17, m = 8.3M B ). Except in the case of 
GJ 3634b, we can note that their masses are large and their radius would also be large in the 
scenario of solid-like planet given in Valencia et al. (2007a). However in all cases the prolateness 
are e <C 0.01 and dynamics of rotation is regular, like several other cases shown previously, in 
spite of their eccentric orbit. 

The case of GJ 3634b is different due to its proximity to GJ 3634 and relatively small 
mass. In this case, due to the similarities of the values of the parameters of the star-planet 
system with the system GJ 876d, we can use the results of the Section 3.2.1 to infer that the 
dynamics of rotation of is regular in spite of a probably non-negligible value of prolateness. 



3.4 Other cases: GJ 1214b and GJ 581g 

GJ 1214b has a mass within the interval m < lOMg, but as shown in Charbonneau et al. 
(2009) and Valencia (2011), due its large radius, it may have a substantial gas layer and 
cannot be considered a solid-like planet. Our model is not adequate to this kind of structure 
but we speculate on the dynamics of rotation of GJ 1214b on basis of our model. And we also 
use this planet to give an application of the Chirikov's overlap criterium. 

Figure 8 shows four plots where we consider our estimative value of e = 0.0186 (Table 
II). Four values of eccentricity, beginning with e = 0.27 (the maximum admitted value by 
Charbonneau et al. 2009), has been used to compute the sections. The results showed in 
Figure 8(a) shows a large region of chaotic motion close to the main resonancef^l. For smaller 
values of eccentricity, the chaotic regions are confined close to the separatrices (e.g. Figure 
8(b)). 



12 Our numerical experiments show that in the cases of strong perturbation, the domain of 
the 1:2 resonance is preserved from the chaotic regions located close to it. Wisdom (1984), in 
their study on the rotation of the satellite Hyperion, discuss this property of the 1:2 resonance. 
In fact, analyses of the averaged torques associated to the 1:1 and 1:2 spin-orbit resonances 
show that they are the strongest ones among all low-order resonances (e.g. Goldreich and Peale 
1966). 
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(de/dt)/n 



e=0.001 



55 Cnc e: (B-A)/C~0.028 
b) 



e=0.01 




Fig. 7 Surfaces of Section in the plane (0/n X 9) of numerical solutions of JTJ. The sections 
have been done at each orbital revolution of the planet 55 Cancri e adopting different values 
of the orbital eccentricity which are indicated at the top of the plots. Each orbit has been 
integrated for 1140 orbital periods of the planet. The current (published) value of eccentricity 
is considered in (c) while the other values are fictitious. -2^1 ~ e = 0.028 are estimated 
adopting R = I.YIKe for the radius of the planet, (e) Dynamical map (see Section 3.2.2) in a 
grid of 22,801 points corresponding to (c). 
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GJ 1214b: (B-A)/C~0.0186 
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Fig. 8 Surfaces of Section in the plane (6/n X 0) of numerical solutions of |T||. The sections 
have been done at each orbital revolution of the planet GJ 1214b, for four different values of 
the orbital eccentricity within the interval e < 0.27 (Charbonneau et al. 2009). Each orbit has 
been integrated for 1140 orbital periods of the planet. 



The other plots given in Figures 8(c,d) correspond to smaller values of the eccentricity and 
show that the dynamics is regular and the domain of the 1:1 resonance is the largest one in 
the phase space. 

We can use the Chirikov's overlap criterium to estimate local chaos around the main 
resonances in terms of the main parameters of the problem: B q A and eccentricity. In fact, 

the separatrices of 1:1 and 3:2 resonances overlap for all values of B q A larger than a critical 
quantity which is given by the following expression: 



(B - A 



) - 


<-( 




' critical 


3 V2 H 





2 



(Wisdom et al. 1984). 

Figure 9 shows the plot of the critical ( B ~c ) • i as a f unc ti°ri of eccentricity. For 
e = 0.27, the critical value is 0.022. We can see that this value is very close to e = 0.0186, 
which has been calculated for GJ 1214b. Thus, in this case the overlap occurs, what is in well 
agreement with the result obtained in the surfaces of section (Figure 8(a)). For smaller values 
of eccentricity (e = 0.1 for instance), B ^ A would be larger than 0.0186 for overlap to occur 
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according to the Equation]?] and for this reason the overlap of the 1:1 and 3:2 resonances does 
not occur in the surface of section (Figures 8(b)). 

We can make an analysis similar to this given above for all other cases shown in this paper. 
For instance, in Figure 9 we indicate the value of e = 0.028 calculated for 55 Cnc e, and we 
have that, for e = 0.15 the overlap occurs. Figures 9(c,d) show good agreement between the 
surfaces of section and the Chirikov's overlap criterium: for e < 0.1 the separatrices of the 1:1 
and 3:2 do not overlap (Figures 7(b,c)) but for e = 0.1 (Figure 7(d)), a valor very close to the 
critical one, the overlap occurs. 

Figure 9 shows that for large values of B ^ A ~ 10" 1 the overlap may occur for very small 
value of eccentricity ~< 0.01. This is shown in Figure 3(c) in the case of KOI-55c. 



GJ 581 g 

We give here some comments about the dynamics of rotation of GJ 581 g. This planet, 
although not yet surely confirmed (see Tadeu dos Santos et al. 2012), would have a very small 
m ~ 3.1Mg, which would imply in a relatively large value of e. However, due to its relatively 
large distance from GJ 581 (r ~ 0.14 AU), and the slow-mass star, its prolateness is negligible, 
and its rotation is regular within the domain of our model. 



4 Conclusions and Discussion 



We have explored the dynamics of rotation of close-in exoplanets with a model often applied 
in studies of physical libration of natural satellites and Mercury (Goldreich and Peale 1966, 
Henrard 1985, Bills et al. 2009). This model is suitable to quantify the main effects of the 
gravitational torque of the star on the rotation of the planet when the obliquity of the axis of 
rotation of planet is considered fixed (i.e., we neglect the motion in attitude). 

The model is valid in the hypothesis of a rigid planet and therefore we neglect in our 
calculations the effects of tidal torques of the central star on the rotation of the planet. Thus, we 
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study the conservative dynamics of rotation of several possible Earth-like cxoplanets. This type 
of planet, which mass is m < 10-Mg, may admit a permanent solid structure depending also on 
the value of the radius of the planet (e.g. Valencia et al. 2007a, 2011). Our methodology was 
based in previous works like Wisdom et al. (1984, 1987, 2004), where the dynamics of rotation of 
regular satellites of Saturn and Phobos and Deimos are numerically explored through analysis 
of surfaces of section. 

According to the conservative model, the effects of the torque of the star on the rotation of 
the planet depends basically on two parameters: the eccentricity of the orbit, and the existence 
of a permanent equatorial bulge (the prolateness, e) of an ellipsoid with semi-axes a > b > c. A 
non-null value of the prolateness implies that B q A 7^ 0, where A, B, C are the main moments 
of inertia which, in the case of exoplanets, arc completely unknown quantities. We show in 
this work that B ^ A = e, a relation valid at first order in e when b = c. On the other hand e 
can be written as function of the masses of the system, the star-planet distance and the radius 
of the planet, and when these parameters arc known, we can estimate an order of magnitude 
of e. In general the transiting planets have a determination of the radius. When the planet is 
not a transiting one we show here how to obtain a significant range of the planet prolateness 
in terms of the several parameters of the star-planet system. 

We discuss the cases of Super-Earths candidates immersed in multi-planetary systems and 
those which are single planets. We also consider planets in circular and eccentric orbits. In the 
case of circular orbits the dynamics of rotation of exoplanets is trivial, i.e., it is regular and the 
synchronous state governs the phase space with a domain dictated only by the magnitude of 
the prolateness of the equator and the mass of the central star. The most interesting cases are 
closc-in exoplanets with eccentric orbits, since in this case the dynamics of physical libration 
may deviate significantly from the pendulum analogy. Large-scale chaotic motion of rotation for 
initial conditions close to the separatrix of the main resonances appears due to a combination 
of large prolateness and high eccentricity. 

Our main results on the dynamics of rotation of some cases studied here are summarized 
below: 

• KOI-55b,c, closc-in planets with very small mass (m < lMg) have large prolateness 
0.1 < e < 0.2 and, even for small eccentricity (e ~ 0.007), the separatrices of the 1:1 and 3:2 
resonance overlap leading to the possibility of complicated dynamics of rotation of the planet 
around the synchronous motion. 

• The rotation of planets like Kepler-lOb CoRoT-7b, KOI-961c can suffer strong effects of 
the star due to their relatively large values of their prolateness, which are given in the interval: 
0.01 < e < 0.02. However, their orbits are probably circular and the phase spaces of rotation 
are dominated by the circulatory regimes and libration around the synchronous island, which 
domain is large and is probably the loci of the current state of the rotations of this kind of 
planet. 

• Planets with prolateness e <C 0.01 and quasi-circular orbits have the phase space dom- 
inated by circulation regimes of motion, and the width of the synchronous resonance is very 
small in the phase space. As examples we can cite CoRoT-7c and KOI-961b. 

• Closc-in planets in eccentric orbits and orbital period P < 7.5 days admit complicated 
dynamics depending on the value of the radius (R) and the mass (m). In fact, the larger is R, 
the larger is the prolateness; the larger is the m, the smaller is the prolateness. For instance, 
GJ 876d can have large prolateness e > 0.01 if R ~ 3Re- It is improbably that 61 Vir b has 
non-negligible prolateness in the case of R 2> 3Re- Kepler-9d, Gl 581e, HD 215487b have 
small prolateness e <C 0.01 and the rotation phase spaces are regular. For eccentricities larger 
than 0.1, we have observed that in some cases (e.g. HD 215487b, Kepler-9d) the domains of 
the 3:2 and 1:1 resonances have very similar widths in the phase space. 

• 55 Cnc-e admits a large chaotic regions close the synchronous rotation adopting its 
current parameters (e = 0.057, m ~ 8.58Mb, R = 2.17-Re, e ~ 0.028). The implications of the 
chaos on the evolution through spin-orbit resonance are being investigated in another work 
(Rodriguez et al. 2012). 

All simulations show that the dynamics of rotation of close-in eccentric Earth-like planets 
may have a complex behavior close to the regions of separatrices in the phase space depending 
on the prolateness. In a scenario including the action of a dissipative force which acts to spin 
down the primordial rotation, one should understand how planetary rotation of all bodies stud- 
ied here evolves to synchronism traversing a chaotic region of the phase space (e.g. Rodriguez 
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et al. 20120 Even in those cases with almost circular orbits, the rotation can be driven to a 
complex motion when — ^ — is large (e.g. KOI-55b,c). 

• The dynamics of planets with eccentric orbits and (P > 7.5 days) have also been studied 
here (e.g. Gl 581c, HD 181433b, HD 215487c). We have shown that all of them have negligible 
prolateness in the case we suppose a solid structure for them. Therefore they probably have 
regular motion of rotation around the synchronous or other main spin-orbit resonances. 

• The prolateness of several single Earth-like candidates planets considered here (e.g. GJ 
176b, GJ 674b, HD 97658b) may be negligible due to both their large values of mass and 
distance from the star. Therefore their dynamics of rotation are regular in spite of eccentric 
orbits (GJ 674b, HD 97658b). On the other hand the planet GJ 3634b has probably a non- 
negligible prolateness e ~ 10~ 3 . With this value as reference, Rodriguez et al. (2012) show 
simulations where this planet may be currently captured into the 2:1 spin-orbit resonance. 



Discussion on Tidal Evolution 

• Due to tidal effects it is expected that the orbit of a close-in planet be circular. In Appendix 
A we calculate the timescale for circularization of a planet with elements similar to CoRoT-7b. 
We confirm previous results which show that state that planets with orbital period P > 8 
days may admit a non-circular orbit since the time scale for circularization is very large: > 10 
billion years. Since the value of orbital eccentricity of some planets are still being improved 
(see Tadeu dos Santos et al. 2012) we consider in our calculations several values of eccentricity 
in addition to the ones given in literature. 

• We can also apply our calculations given in Appendix A for several single planets dis- 
cussed in Section 3.4. The orbital periods of GJ 3634b, GJ 674b, HD 7924b and HD 45184b 
are smaller than 6 days and thus it is expected a circular orbit in view of the tidal effect. On 
the other hand, it is not be expected a null eccentricity for HD 125595b since P ~ 9.4 days. 

• Although we have neglected the effects of tidal torques in this work, they would be 
important in the case of a non-rigid planet located very close to the central star. In Appendix 
B we study the case of the planet CoRoT-7b utilizing a model of tides where the phase lags 
depend linearly with the corresponding tidal frequency. We calculate the magnitude of the 
tidal torque and compare with the gravitational torque due to the equatorial prolateness in 
order to evaluate the stability of two spin-orbit resonances, the synchronous and 3:2 spin-orbit 
resonances. The 1:1 resonance is stable adopting our estimative of e = 0.0092 for CoRoT-7b. 
The reader is referred to Rodriguez et al. (2012) for a detailed study of the coupled spin-orbit 
evolution of short-period rocky planets. 
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Appendix A: Orbital circularization due to tidal effect 

We consider a slow-rotating star and a close-in planet interacting pair. The aim is to analyze 
the timescale for orbital circularization due to tidal interaction. We refer to Ferraz-Mello et al. 
(2008) and Rodriguez and Ferraz-Mello (2010) for assumptions, definitions of quantities and 
further details. 

The averaged variation of the eccentricity due to the combined effects of stellar and plan- 
etary tides is given by 




nea 5 (18s + 7p), 



(8) 



where 



4 Qo rn 




(9) 
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are two parameters which stand for stellar and planetary tides, respectively. The symbol 
refers to the star, k 2 is the second degree Love number, Q is the dissipation function or quality 
factor. It can be shown that the orbital circularization can be accounted by planetary tides 
only. Indeed, § is proportional to (m/mo)(l/Qo) which becomes a small quantity for small 
mass planets and typical stellar Qo values. Thus, the contribution of stellar tides can be safely 
neglected in our analysis (see Rodriguez and Ferraz-Mello 2010). 

The timescale for orbital circularization can be denned by r e = e/|(e)|, or 



1„5 



3n _1 a' 

Writing Equation (1101) as a function of the semi-major axis, we obtain 



(10) 



Aa 13 / 2 , (11) 

ernatively, we can ex 

P in the following way 



where A = 3 f . Alternatively, we can express the result as a function of the orbital period 

7v Gmn 



Te =BP 13/3 , (12) 

/ \13/6 

where B = A I ™° ) • Figure 1101 shows the plot of Equation (1 1 2 1> for a planet with the 

properties of CoRoT-7b, a Super-Earth planet with m = 8.5m E (Ferraz-Mello et al. 2011), 
R = 1.68i?, B , assuming the values Q = 100, Q = 500, Q = 1000, and k 2 = 0.35. 

We clearly see that the circularization timescale decreases for short-period planets, as 
Equation il 1 21) indicates. Moreover, noting that r e <x Q, the orbital circularization would be- 
comes faster in the case of small Q values (i.e., large dissipation). As an example, for P = 4 
days and Q = 100 we have r e ~ 181 Myr. 

It is important to note that the parameter B is linearly dependent on Q/k 2 , which is 
a quantity poorly known for extrasolar planets. Hence, the plot shown in Figure 1101 can be 
strongly modified if others values of the planet dissipation are considered. 

The result (1 1 2[ i can be useful in order to quantify the efficiency of tides to produce orbital 
circularization of close-in planets. Note that, in some cases, r e can be compared with the age 
of the system, indicating that the orbit of the close-in planet should be circularized during the 
planet lifetime. However, for large P, T e can be even larger than the age of a typical planetary 
system, in which case a non-circular orbit should be expected due to tidal interaction. 



Appendix B: Capture in spin-orbit resonance 

The numerical exploration of the planet rotation, which is subject to the gravitational torque of 
the star, has shown different behaviors including the oscillation around spin-orbit resonances. 
When a dissipative effect like the tidal torque is included, the rotation can be captured in a 
resonant motion. The specific capture depends on the eccentricity, and also on Q and e. Hence, 
as e is tidally damped, the capture should become unstable and the rotation can achieves 
another resonant state, which, at the same time, should result in a temporary trapping. When 
the orbit completes the circularization due to the tidal torque, the final evolution result in the 
synchronization between the orbital and rotational periods (i.e., the 1:1 spin-orbit resonance). 
The reader is referred to Goldreich and Peale (1966) and Rodriguez et al. (2012) for further 
details on the spin-orbit evolution of close-in planets. 

Let us first consider the torque due to the prolateness or permanent equatorial deformation 
(i.e., a / b) on a rotating body of mass m and radius R. The maximum torque, averaged over 
an orbital period is given by 

< N >= -^n 2 {B - A)H(p,e) (13) 

(see Goldreich and Peale, 1966 and Equation Q), where H(p,e) are power series in e and 
p = fi/n, with 17 the angular velocity of rotation of the deformed body (Goldreich and Peale 
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Fig. 10 Plot of Equation H21 for CoRoT-7b super-Earth planet. The mass of the host star 
is mo = 0.93mg(7jv 



1966). Equation (1136 assumes that there is a commensurability between Q and n, indicating 
that p= 1,-1/2,1,1/2... 

In addition to the above torque, we also consider the tidal torque driven by the central 
body of mass mg. The averaged tidal torque reads 



< T >= ^-T I 4e o + e 2 (-20e + 49 £l + e 2 (14) 

8a° 

(see Ferraz-Mello et al. 2008), where &2 is the second degree Love number, whereas Ei are 
the phase lags of tidal waves with frequency m. The phase lags account for the internal vis- 
cosity, which introduces a delay between the action of the tidal force and the corresponding 
deformation. 

Several tidal models can be used to fix the dependence between phase lags and frequencies, 
that is, the function £i(vi). We first consider the usually called linear model, where Ei = ViAt, 
where At is known as time lag and is considered constant in the linear model (Mignard, 1979). 
The frequencies associated to the phase lags appearing in Equation (1 1 4 H are uo = 2Q — 2n, 
v\ = 2Q — 3n and = 21? — n (see Ferraz-Mello et al. 2008). Replacing into Eq. J 141) and 
applying the linear model, we obtain 



3k 2 AtGmlR b n „ 
<T>= [2p-2 + (15p-27)e 2 ]. (15) 

The time lag can be related to the most used quantity Q, the dissipation function. Since 
Si ~ Q" 1 , it follows, under assumption of linear model, Qq~ ~ 2nAt(p— 1), Q^ 1 ~ nAt(2p— 3), 
Q^ 1 ~ nAt(2p - 1). We note that singularities in Qi are associated to spin-orbit commensu- 
rability with p = 1, p = 3/2 and p = 1/2. Since we restrict our investigation to the cases of 
1 : 1 and 3 : 2 spin-orbit resonances, we call Q = Q2 in order to avoid (^-singularities. Hence, 
the relationship between Q and At is given by Q = l/[nAt(2p — 1)]. 
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Fig. 11 Critical value of the equatorial ellipticity as a function of orbital eccentricity for two 
spin-orbit resonances. The CoRoT-7+7b system illustrates the example. 



Stationary solutions 

By definition, the stationary solutions of the rotation are those which satisfy 

<AT> + <T>=0. (16) 

Because we are interesting in those solutions for which there exist commensurability between 
spin and orbital revolutions, we can determine a critical value of e (we call it e*) which allows 
a spin-orbit resonance motion be maintained when the planet's rotation is under simultaneous 
action of two torques. Thus, using Equation H13I15I161 we obtain 

1 k 2 m /R\ 3 [2p - 2 + (15p - 27)e 2 ] 

I Q m {a J (2p-l)H(p,e) ' ( ' 

where we have used the third Kepler law and (B — A) ~ Ce = £mR 2 e, where C is the moment 
of inertia about the rotation axis and < £ < 2/5. 

The condition e > e* is usually known as stability condition of the p- resonance (e.g. 
Goldreich, 1966). In fact, the stability condition requires that < T > not exceed the maximum 
restoring torque < N > and, for that reason, e* should be considered as a critical value. 

We note that the case p = 1 is in agreement with the result found in Ferraz-Mello et al. 
(2008) for the synchronous motiorP^I. 

Figure [TT] shows the variation of e* with the orbital eccentricity, taking the planet CoRoT- 
7b as an example. We adopt k 2 = 0.35, £ = 0.35 and Q = 100. The cases of 1 : 1 (synchronous 
rotation) and 3 : 2 spin-orbit resonances are shown. To second order in eccentricity, we have 
H (1, e) = 1 — 5e 2 /2 and H(3/2,e) = 7e/2. The dashed horizontal line indicates the value of 
e used in the simulations for CoRoT-7b (e = 0.00992; Figures 2(a,b)). We note that e > e* in 
both cases, indicating that the resonant motion should be stable for the range of considered 
eccentricity. However, e* can reach high values for very small e as Equation H17II indicates. On 
the other hand, we have seen in the numerical simulations (no tides case, Figure 2(b)) that 
the domain of the 3 : 2 spin-orbit resonance is very small for almost-circular orbits. 



14 For comparison it is necessary to consider the relationship between J22 (the equatorial 
ellipticity) and B — A, that is, B - A = 4J 2 2mR 2 (see Beutler 2005). 



